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Eigenconfigurations of Tensors 

Hirotachi Abo, Anna Seigal, and Bernd Sturmfels 

Abstract. Square matrices represent linear self-maps of vector spaces, and 
their eigenpoints are the fixed points of the induced map on projective space. 
Likewise, polynomial self-maps of a projective space are represented by tensors. 
We study the configuration of fixed points of a tensor or symmetric tensor. 


1. Introduction 

Square matrices A with entries in a field K represent linear maps of vector 
spaces, say AT" —)• A"", and hence linear maps ip : p"-i of projective 

spaces over K. If A is nonsingular then ■)/) is well-defined everywhere, and the 
eigenvectors of A correspond to the fixed points of ijj. The eigenconfiguration of A 
consists of n points in P"“^, provided A is generic and K is algebraically closed. 
Conversely, every spanning configuration of n points in P"“^ arises as the eigen¬ 
configuration of an n X n-matrix A. However, for special matrices A, we obtain 
multiplicities and eigenspaces of higher dimensions [SI. Moreover, if A" = R and 
A is symmetric then its complex eigenconfiguration consists of real points only. 

This paper concerns the extension from linear to non-linear maps. Their fixed 
points are the eigenvectors of tensors. The spectral theory of tensors was pioneered 
by Lim [Limj and Qi |Qi| . It is now a much-studied topic in applied mathematics. 
For instance, consider a quadratic map ip : P"“i P"“i, with coordinates 

n n 

(1.1) 1pi{xi, . . . ,Xn) = EE QijkXjXk for i = 1,..., n. 

i=i fe=i 

One organizes the coefficients of ip into a tensor A = (oyfe) of format n x n x n. 

In what follows, we assume that A = (aqia-.-id) is a d-dimensional tensor of 
format nxnx ■ ■ ■ xn. The entries lie in an algebraically closed field AT of 

characteristic zero, usually the complex numbers AT = C. Such a tensor A € 
defines polynomial maps A" —>■ A" and p”-i just as in the formula (HI: 

n n n 

ipi{xi,...,Xn) = E E ■ ■ ■ E for i = 1,..., n. 
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Thus each of the n coordinates of '0 is a homogeneous polynomial 'ipi of degree d—1 
in X = (xi,X 2 , ■. ■ ,Xn)- The eigenvectors of A are the solutions of the constraint 


«x) ::: «x)) ^ 

The eigenconfiguration is the variety defined by the 2 x 2-minors of this matrix. 
For a special tensor A, the ideal defined by (HID may not be radical, and in that 
case we can study its eigenscheme. Recent work in develops this for d = 2. 

We note that every n-tuple (01,..., 0„) of homogeneous polynomials of degree 
d — 1 in n variables can be represented by some tensor A as above. This repre¬ 
sentation is not unique unless we require that A is symmetric in the last d—1 
indices. Our maps ip : — ->■ p"-i are arbitrary polynomial dynamical system 

on projective space, in the sense of |FS| . Thus the study of eigenconfigurations of 
tensors is equivalent to the study of fixed-point configurations of polynomial maps. 
Of most interest to us are symmetric tensors A, i.e. tensors whose entries 
8-re invariant under permuting the d indices. These are in bijection with 
homogeneous polynomials (p = 8nd we take ipj = d(p/dxj. 

The eigenvectors of a symmetric tensor correspond to fixed points of the gradient 
map V(p : P"“i P”“i, and our object of study is the variety in P"“i defined by 


/-| Q\ f 

^ \d(p/dxi d(p/dx2 d(p/dxn 

This paper uses the term eigenpoint instead of eigenvector to stress that we 
work in P"“^. In our definition of eigenpoints we include the common zeros of 
01,... ,0„. These are the points where the map P”“^ —p"-i is undefined. For 
a symmetric tensor 0, they are the singular points of the hypersurface {0 = 0} in 
P"“^. At those points the gradient V0 vanishes so condition (11.31) holds. 



Example 1.1. Let n = d = 3 and 0 = xyz. The corresponding symmetric 
3x3x3 tensor A has six nonzero entries 1/6 and the other 21 entries are 0. Here 
V0 : P^ ---> P^, (a: : 2 / : z) ^ {yz : xz : xy) is the classical Cremona transformation. 
This map has four fixed points, namely (1 : 1 : 1), (1 ; 1 : —1), (1 : —1 : 1) 
and (—1 : 1 : 1). Also, the cubic curve {0 = 0} has the three singular points 
(1 ; 0 : 0), (0 : 1 ; 0), (0 : 0 ; 1). In total, the tensor A has seven eigenpoints in P^. 


This paper is organized as follows. In Section 2 we count the number of eigen¬ 
points, and we explore eigenconfigurations of Fermat polynomials, plane arrange¬ 
ments, and binary forms. Section 3 generalizes the fact that the left eigenvectors 
and right eigenvectors of a square matrix are distinct but compatible. We explore 
this compatibility for the d eigenconfigurations of a d-dimensional tensor with n = 2. 

Section 4 concerns the eigendiscriminant of the polynomial system and 
its variant in (lO) . This is the irreducible polynomial in the n‘^ unknowns 
which vanishes when two eigenpoints come together. We give a formula for its 
degree in terms of n, d and i. Section 5 takes first steps towards characterizing 
eigenconfigurations among finite subsets of P”“^, starting with the case n = d = 3. 

In Section 6 we focus on real tensors and their dynamics on real projective space 
Pg~^. We examine whether all complex eigenpoints can be real, and we use line 
arrangements to give an affirmative answer for n = 3. The paper concludes with a 
brief discussion of attractors for the dynamical systems 0 : --■» These 

are also known as the robust eigenvectors of the tensor power method 
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2. The count and first examples 


In this section we assume that the given tensor A is generic, meaning that it lies 
in a certain dense open subset in the space of all nx ■ ■ ■ xn-tensors. This set 

will be characterized in Section 4 as the nonvanishing locus of the eigendiscriminant. 


Theorem 2.1. The number of solutions in P” ^ of the system equals 


( 2 . 1 ) 


(d- 1)" - 1 
d-2 




i^O 

The same count holds for eigenconfigurations of symmetric tensors, given by El). 
In the matrix case (d = 2) we use the formula on the right, which evaluates to n. 


This result appeared in the tensor literature in [CSlIOO], but it had already 
been known in complex dynamics due to Fornaess and Sibony [FSl Corollary 3.2]. 
We shall present two proofs of Theorem l2.ll cast in a slightly more general context. 

For certain applications (e.g. in spectral hypergraph theory |LQY| ), it makes 
sense to focus on positive real numbers and to take the root after each iteration of 
the dynamical system if. This leads to the following generalization of our equations: 


( 2 . 2 ) 


/ x( X 

V'0l(x) V’2(x) 



< 1 . 


We refer to the solutions as the eigenpoints of the given tensor A. For i = 1, this 
is the definition in the Introduction. In the nomenclature devised by Qi [CQZKog, 
one obtains E-eigenvectors for £ = 1 and Z-eigenvectors for £ = d—1. The subvariety 
of P"“i defined by (12.21) is called the eigenconfiguration of the tensor A. 


Theorem 2.2. 
(2.3) 

distinct points in P 


The eigenconfiguration of a generic tensor A consists of 

n—1 


(d- l)”" - 


d-1 - £ 


i=0 


n — 1 


. If £ = d — 1 then the formula on the right is to he used. 


Proof. Consider the 2xn-matrix in (|2.2I) . Its rows are filled with homogeneous 
polynomials in S' = K[xi, ..., Xn] of degrees £ and m respectively, where the tpi are 
generic. Requiring this matrix to have rank < 1 defines a subscheme of P”“^. By the 
Thom-Porteous-Giambelli formula [m § 14.4], this scheme is zero-dimensional, 
and its length is given by the complete homogeneous symmetric polynomial of 
degree n — 1 in the row degrees, £ and m. This is precisely (ESI) if we set m = d—1. 

Another approach, which also shows that the scheme is reduced, is to use vector 
bundle techniques. Consider the 2 x n-matrix as a graded S'-module homomorphism 
from S{—£) © S{—m) to S'®". The quotient module Q of S®" by the submodule 
generated by the first row {x\,... ,xlf) is projective. In other words, the sheafifi- 
cation Q of Q is locally free. The scheme associated with the 2 x n-matrix can 
therefore be thought of as the zero scheme of a generic global section of Q{m). 
Since Q{m) is globally generated, the scheme is reduced |Ein[ Lemma 2.5]. □ 

Here is a brief remark about eigenvalues. If x £ df" is an £^^ eigenvector of A 
then there exists a scalar \ G K such that = Xxf for all i. We call (x. A) an 

eigenpair. If this holds then is also an eigenpair for all v £ A'\{0}. 

Such equivalence classes of eigenpairs correspond to the £^^ eigenpoints in P”“^. 
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The case i = d—1 is special because every eigenpoint has an associated eigenvalue. 
If ^ ^ d — 1 then eigenpoints make sense but eigenvalues are less meaningful. 


Proof of Theorem 12.11 The first statement is the £= I case of Theo¬ 
rem 12.21 For the second assertion, it suffices to exhibit one symmetric tensor ip 
that has the correct number of eigenpoints. We do this for the Fermat polynomial 

(2.4) ~ Xi F X 2 F ■ ■ ■ F x'^. 

According to ( 121 ), the eigenconfiguration of (p is the variety in P" ^ defined by 


(2.5) ■■■ 2i) < 1. 

We follow [Rob] in characterizing all solutions x in to the binomial equations 

XiXj{x^~‘^ — Xj~'^) = 0 for 1 < i < j < n. 

For any non-empty subset I C {1,2,. ..,n}, there are (d — 2)1'^!“^ solutions x with 
supp(x) = {i \ Xi ^ 0} equal to I. Indeed, we may assume Xi = \ for the smallest 
index i in /, and the other values are arbitrary (d — 2)'^'^ roots of unity. In total. 


^(d-2)l^l-i 




(d-2-h 1)" - 1 

dPl 


This equals m- Here we assume d > 3. The familiar matrix case is d = 2. □ 


Example 2.3. Let d = 4. For each /, there are 21-^1 ^ eigenpoints, with Xi = ±1 
ior i G I and Xj = 0 for j ^ I. The total number of eigenpoints in P"“i is (3" —1)/2. 

We note that the argument in the proof of Theorem l2.1l does not work for i >2. 
For instance, if .^ = d — 1 then every point in P"“^ is an eigenpoint of the Fermat 
polynomial. At present we do not know an analogue to that polynomial for i > 2. 


Problem 2.4. Given anyi,d andn, exhibit explicit polynomials ^(x) of degree 
d in n variables such that i2.2\) has \2. d)) distinct isolated solutions in P"“^. 


We are looking for solutions with interesting combinatorial structure. In Sec¬ 
tion 6 we shall examine the case when (/>(x) factors into linear factors, and we shall 
see how the geometry of hyperplane arrangements can be used to derive an answer. 
A first instance was the Cremona map in Example ll.il Here is a second example. 

Example 2.5. For n = 4 the count of the eigenpoints in (12.11) gives d^—2d^-|-2d. 
We now fix d = 5, so this number equals 85. Consider the special symmetric tensor 

(p{x) = XiX2X3X4{Xi F X2 F X3 F X4:). 

The surface defined by <p consists of five planes in P^. These intersect pairwise in 
ten lines. Each point on such a line is an eigenpoint because it is singular on the 
surface. Furthermore, there are 15 isolated eigenpoints; these have real coordinates: 

(2 : 2 : -1 : -1), (2 : -1 : 2 : -1), (2 : -1 : -1 : 2), (-1:2:2: -1), 

(2.6) (-1 : 2 : -1 : 2), (-1 : -1 : 2 : 2), (1:1:1: 1), (i(5 ± ^TS) : 1 : 1 : l), 

(1 : i(5±yT3) : 1 : 1), (1 : 1 : i(5±yT3) : l), (l : 1 : 1 : i(5±yT3)). 

The five planes divide P^ into 15 regions. Each region contains one point in (12.61) . 

Now, take a generic quintic (p'{x.) in R[a:i,a;2,cca,0:4], and consider the eigen¬ 
configuration of (ppx.) Fe(p'{x.). This consists of 85 points in P^. These are algebraic 
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functions of e. For e > 0 small, we find 15 real eigenpoints near (IT^ . The other 70 
eigenpoints arise from the 10 lines. How many are real depends on the choice of (j)'. 

The situation is easier for n = 2, when the tensor A has format 2x2x • • • x2. 
It determines two binary forms ipi and tjj 2 - The eigenpoints of A are defined by 

(2.7) y'tl}i{x,y) - X ■ ip2{x,y) = 0. 

This is a binary form of degree d, so it has d zeros in as predicted by (EU. 
Conversely, every binary form of degree d can be written as yipi—xil} 2 - This implies: 

Remark 2.6. Every set of d points in is the eigenconfiguration of a tensor. 


The discussion is more interesting when we restrict ourselves to symmetric ten¬ 
sors. These correspond to binary forms (p{x^ y) and their eigenpoints are defined by 


d(j) dcj) 


0 . 


The matrix case (d = 2) shows that r{emark [2.6l cannot hold as stated for symmetric 

tensors. Indeed, A A = and </> = ax'^ + 2bxy + cy'^ then = 

—bx^ -|- (a — c)xy + by^. This confirms the familiar facts that the two eigenpoints 
(ui : ui) and {u 2 '■ V 2 ) are real when a,6,c G R and they satisfy U 1 U 2 + V 1 V 2 = 0. 
The following result generalizes the second fact from symmetric matrices to tensors. 


Theorem 2.7. A set of d points {ui : Vi) in P^ is the eigenconfiguration of a 
symmetric tensor if and only if either d is odd, or d is even and the operator 

/ q2 q2 \d/2 d- 

2 + 'g~2 j annihilates the corresponding binary form — Uiy). 


Proof. The only-if direction follows from the observation that the Laplace 
operator ^ commutes with the vector field y-^ — x-^. Hence, for any cj) of 

degree d, we obtain zero when ^2 + ^ gets applied d/2 times to 

For the if direction, we examine the (d -|- 1) x (d -|- l)-matrix that represents 
the endomorphism 4> y^ space of binary forms of degree d. This 

matrix is invertible when d is odd, and its kernel is one-dimensional when d is even. 
Hence the map is surjective when d is odd, and it maps onto a hyperplane when d 
is even. The only-if part shows that this hyperplane equals ( -I- )^- D 

After completion of our manuscript we learned that Theorem l2.7l was also found 
independently by Mauro Maccioni, as part of his PhD dissertation at Firenze, Italy. 


Example 2.8 (d = 4). Eour points (mi:ui), {u 2 .V 2 ), (us'.vs), {u 4 :v 4 ) on the line 
P^ arise as the eigenconfiguration of a symmetric 2 x 2 x 2 x 2-tensor if and only if 

3U1U2U3U4 + U1U2V3V4 -I- UlU3Z;2U4 -I- U1U4U2U3 -I- • • • -I- U3,U4VlV2 + 3uiU2U3U4 = 0 . 


This equation generalizes the orthogonality of the two eigenvectors of a symmetric 
2 X 2-matrix. For instance, the columns ofC/=fQ | M represent the eigen¬ 


configuration of a symmetric 2 x 2 x 2 x 2-tensor, but this does not hold for 


); 


u. 
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Example 12.81 underscores the fact that the constraints on eigenconfigurations 
of symmetric tensors A are not invariant under projective transformations. They 
are only invariant under the orthogonal group 0(n), like the Laplace operator in 
Theorem 12.71 By contrast, the constraints on eigenconfigurations of general (non- 
symmetric) tensors, such as Theorem 15.11 will be properties of projective geometry. 

We are familiar with this issue from comparing the eigenconfigurations of real 
symmetric matrices with those of general square matrices. These are respectively 
the 0(n)-orbit and the GL(n)-orbit of the standard coordinate basis. 


3. Compatibility of eigenconfigurations 


When defining the eigenvectors of a tensor A, the symmetry was broken by 
fixing the first index and summing over the last d — 1 indices. There is nothing 
special about the first index. For any k G {1,..., d} we can regard A as the self-map 

^[fc] . pn-l pn-1 


whose coordinate is the following homogeneous polynomial of degree d — 1 in x: 


n n 


= E---E E---E 


(l'i^ ... 


+ ' ^jk-l^jk + 1 ' ' ' '^3d' 


• Xn 




Let Eigf^^(^) denote the subvariety of P” ^ consisting of the fixed points of 
For a generic tensor A, this is a finite set of points in P”“i of cardinality 

^ ford >3. 

This raises the following question: Suppose we are given d configurations, each 
consisting of D points in P"“^, and known to be the eigenconhguration of some 
tensor. Under what condition do they come from the same tensor A1 

We begin to address this question by considering the case of matrices (d = 2), 
where D = n. Our question is as follows: given an n x n-matrix A, what is the 
relationship between the left eigenvectors and the right eigenvectors of A1 


Proposition 3.1. Let {vi,V 2 , ..., v„} and {wi,W 2 ,... ,w„} be two spanning 
subsets o/P"“^. These arise as the left and right eigenconfigurations of some nxn- 
matrix A if and only if, up to relabeling, the dot products of vectors corresponding 
to Wi and Vj are zero whenever i j ■ 

Proof. Let U be a square matrix whose columns are the eigenvectors of A. 
Then the columns of form a basis of eigenvectors for . □ 


The condition in Proposition 13.11 defines an irreducible variety, denoted EC „^2 
and called the eigencompatibility variety for n x n-matrices. It lives naturally in 
the space of pairs of unordered configurations of n points in P”“^. In symbols, 

(3.1) EC „.2 C Sym„(P"-i) x Sym„(P"-i). 

It has middle dimension n{n — 1), and it maps birationally onto either factor. We 
may identify Sym„(P”“^) with the Chow variety of products of n linear forms in 
n variables. Here, each configuration {vi, V 2 ,..., v„} is represented by nr=i(''^*'x)- 
The coefficients of this homogeneous polynomial serve as coordinates on Sym„(P"“^) 
It would be worthwhile to express Proposition 13.II in these coordinates. 
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Example 3.2 (n = 2). The eigencompatibility variety EC 2,2 for 2 x 2-matrices 
is a surface in (Sym 2 (P^))^. This ambient space equals (P^)^, by representing a 
pair of unlabeled points on the line P^ with the binary quadric that defines it. To be 
precise, a point ((uo:mi:u 2 )j (■po:^^i:'P 2 )) in (P^)^ is identified with the binary forms 

f{s,t) = uos'^ + uist + U2t^ and g{s,t) = vqs'^ + vist + V2t'^ ■ 

We want the zeros of f{s,t) and g{s,t) to be the right and left eigenconfigurations 
of the same 2 x 2-matrix. Proposition 13.11 tells us that this is equivalent to 

f{s,t) = X{as + bt){cs + dt) and gis,t) = g,{bs — at){ds — ct). 


By eliminating the parameters a, b, c, d, A, and /x, we find that the surface EC 2,2 is 
essentially the diagonal in (P^)^. It is defined by the determinantal condition 


(3.2) 


rank 


Uq 

Ui 

U2 

V2 

-Vi 

Vo 


< 1 . 


Our aim in this section is to generalize this implicit representation of EC 2 , 2 - 


Let EC„_(i denote the eigencompatibility variety of d-dimensional tensors of 
format nxnx ••• xn. This is defined as follows. Every generic tensor A has d 
eigenconfigurations. The eigenconfiguration with index k of the tensor A is the fixed 
locus of the map Each configuration is a set of unlabeled D points in P”“^, 
which we regard as a point in Sym£i(P"“^). The d-tuples of eigenconfigurations, 
one for each index fc, parametrize 

(3.3) EC„.d C (Sym^(P'^-i))''. 

Thus EC„_(i is the closure of the locus of d-tuples of eigenconfigurations of tensors. 

Already the case of binary tensors (n = 2) is quite interesting. We shall sum¬ 
marize what we know about this. Let A be a tensor of format 2 x 2 x • • • x 2, with 
d factors. Each of its d eigenconfigurations consists oi D = d points on the line P^. 
The symmetric power Sym^ (P^) is identified with the P'^ of binary forms of degree 
d. The zeros of such a binary form is an unlabeled configuration of d points in P^. 
Thus, the eigencompatibility variety for binary tensors is a subvariety 

EC2,d C (P^)'^. 

The case d = 2 was described in Example 13.21 Here are the next few cases. 


Example 3.3 (d = 3). Points in (P^)^ are triples of binary cubics 


f{s,t) = UqS^ + Uis'^t + U 2 St'^ + , 

g{s,t) = VqS^ + Vis'^t + V 2 St‘^ + 

h{s,t) = WqS^ + Wis'^t + W 2 St'^ + w^t^, 


where two binary cubics are identified if they differ by a scalar multiple. The three 
eigenconfigurations of a 2x2x2-tensor A = (a^fc) are dehned by the binary cubics 

/(s,t) = A • (a2iis^ — (am — 0212 — a22i)s^t + (0222 — 0112 — ai2i)st^ — 0122^^) 

g(s, t) = fx ■ (ai2is^ — (am — 0122 — a22i)s^t + (0222 — aii2 — a2ii)st^ — a2i2t^) 

h{s,t) = V • (aii2S^ — (am — ai22 — a2i2)s^t -I- (a222 — 0121 — a2ii)st^ — 0221^^) 


Our task is to eliminate the 11 parameters a^-fc and A,/i, p from these formulas. 
Geometrically, our variety EC 2,3 is represented as the image of a rational map 


(3.4) 


(p3)3^ A^(/,g,d). 


.7 
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This is linear in the coefRcients aijk of A and maps the tensor to a triple of binary 
forms. To characterize the image of (13.4|) . in Theorem 13.61 we introduce the matrix 

( Ml - M3 Ml — M3 Mo - M2 Mq - M2 
Ml - M3 0 Mo - M2 0 

0 Wi — 0 Wq — W2 

Let I be the ideal generated by the 3 x 3-minors of E 3 . Its zero set has the 
eigencompatibility variety EC 2,3 as an irreducible component. There are also three 
extraneous irreducible components, given by the rows of the matrix: 

h = {uo - U2,ui - M3), I2 = (vo - M2,mi - M3), and Is = (wg - W 2 ,wi - W3). 



The homogeneous prime ideal of EC 2,3 is found to be the ideal quotient 


(3.6) (/ : / 1 / 2 / 3 ) 


2 X 2-minors of 


/mo - M 2 
V^i - M 3 


Mo - M 2 
Ml - M 3 


Wo - W2 
Wi - W3 


We conclude that the eigencompatibility variety EC 2,3 has codimension 2 in (P^)^. 


Example 3.4 (d = 4). Points in (P"^)^ are quadruples of binary quartics 

Mos"' -I- UiS^t -I- U2S^t^ + Ugst^ + U4t'^, 

Vqs'^ + ViS^t -|- V2S^t^ + V3St^ + M4t^, 

WqS^ + WiS^t + W2S^t^ + W3St^ + W4IA, 

XoS^ + XiS^t + X2s'^t^ + X3St^ + Xit^. 


One can represent the homogeneous ideal of the eigencompatibility variety EC 2,4 in 
a similar way to Example 13.31 Let I be the ideal generated by the 4 x 4-minors of 


UI—U3 

U1—U3 

U1—U3 

U2—U0+U4, 

M2—Mo —M4 

2iio—ii 2 H" 2 ii 4 

3 mo— 2 m 2 + 3 m 4 

CO 

1 

0 

0 

V0—V2+V4, 

0 

V2 

M2 

0 

W3—WI 

0 

0 

IMo-MI2-I-MI4 

W2 

W2 

0 

0 

X3—X1 

0 

0 

Xo + X 4 

X2 


Let lij be the ideal generated by the 2 x 2-minors of the 2 x 7-submatrix consisting 
of the and rows in (13.711 . The homogeneous prime ideal of EC 2,4 C (P^)^ is 
obtained as the ideal quotient (/ : I 12 I 13 I 14 I 23 I 24 I 34 ) ■ We obtain dim(EC 2 , 4 ) = 12. 


Example 3.5 (d = 5). The eigencompatibility variety EC 2,5 has codimension 
4 in (P®)®, so dim(EC 2 , 5 ) = 21. We represent this variety by the 5 x 8 -matrix 


(3.8) 


— Ml -I- U3 — Ms 
Ml — M3 -I- Ms 
—Ml -I- M3 — Ms 
—Ml -I- M3 — Ms 
Uo — U2 + M4 
0 
0 
0 


Ml — M3 -I- Ms 
0 
0 
0 
0 

Mo — M2 -I- M4 
0 
0 


0 

IMl — Mis -I- Mis 
0 
0 

Wo — VJ2 + MI4 

Mlo — MI2 -I- MI4 

IMo — MI2 -I- MI4 

Mlo — MI2 -I- MI4 


0 

0 

Xl — X3 -I- xs 
0 
0 
0 

xo — X2 -I- X4 
0 


0 

0 

yi-y3 + V5 
0 
0 
0 

yo -y 2 + y4. 


As before, the variety of maximal minors of this 5 x 8 -matrix has multiple com¬ 
ponents. Our variety EC 2,5 is the main component, obtained by taking the ideal 
quotient by determinantal ideals that are given by proper subsets of the rows. 


In what follows we derive a general result for binary tensors. This will explain 
the origin of the matrices (13.5L (13.711 and (13.811 that were used to represent EC 2 ,(i- 
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Fix V = K^. Tensors A live in the space For each k, the map A i—>■ 
factors through the linear map that symmetrizes the factors indexed by [d]\{fc}: 

(3.9) F®'" ^ Sym,_i(F)®y, 

where {ei,..., e^} is a basis for V. Taking the wedge product with (xi,X 2 , • ■ •, Xn) 
defines a further linear map 

Sym^-i (F) (g) F -^ Symrf(F) V 


(3.10) UJ 


UJ 


Write for the composition of (13.101) after (13.91) . Thus (Al) is a vector of length 
(”) whose entries are polynomials of degree d that define the eigenconfiguration with 
index k. For instance, in Example 13.31 / = (}d\A), g = £^^^(A), h = (\'^\A). 

The kernel of consists of all tensors whose eigenconfiguration with index 
k is all of P”“^. We are interested in the space of tensors where this happens 
simultaneously for all indices k: 

d 

(3.11) Kn4 = Pi kernel(£f^]). 

k^l 

The tensors in Kn,d can be regarded as being trivial as far as eigenvectors are 
concerned. For instance, in the classical matrix case (d = 2), we have 

Kn,2 = kernel(£f^^) = kernel(£f^^), 

and this is the 1-dimensional space spanned by the identity matrix. 

In what follows we restrict our attention to binary tensors (n = 2). We regard 

l[k] 

as a linear map ^ P'^. The eigencompatibility variety EC 2 ,d is the 

closure of the image of the map (P'^)'^ given by the tuple 

Let be a column vector of unknowns representing points in the fc**' factor P"^. 


Theorem 3.6. There exists a d'Ke-matrixYid with e = dim{K 2 ^d)+d{d+l) — 2‘^, 
whose entries in the fc**' row are Ij-linear forms in such that EC 2 ,d is an 

irreducible component in the variety defined by the d x d-minors of Ed- Its ideal is 
obtained from those d x d-minors by taking the ideal quotient (or saturation) with 
respect to the maximal minor ideals of proper subsets of the rows of . 


Proof. We shall derive this using the linear algebra method in [A^ §2], We 
express £[k] 

as a (d+1) X 2'^-matrix, and we form the d{d -I- 1) x {2'^ (i)-matrix 


(3.12) 


The left d{d -I- 1) x 2‘^-submatrix has entries in {—1, 0, -1-1} and its kernel is K 2 ^d- 
The rank of that submatrix is r = 2^^ — dim(iF 2 ,d). Using row operations, we can 

.0 u, 

C is an e X d-matrix whose fc*** column has linear entries in the coordinates of . 



uW 

0 ■ 

■ 0 \ 

im 

0 

u[2] . 

0 

yild] 

0 

0 • 

• ) 


transform (13.121) into a matrix 


where A is an r x 2^^ matrix of rank r, and 
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The variety EC 2 ,(i is the set of all points , ul'^1) in such that 

the kernel of (13.1211 contains a vector whose last d coordinates are non-zero. Equiv¬ 
alently, the kernel of C contains a vector whose d coordinates are all non-zero. 

Let Ed be the transpose of C. This is a d x e-matrix whose row has entries 
that are Z-linear in u^^l. By construction, EC 2 ,d is the set of points ..., 
in such that v-E^ = 0 for some v G (if\{0})'^. This completes the proof. □ 


By our matrix representation, the codimension of EC2,(i is at most e — d-|-1, so 
(3.13) dim(EC2,d) > — (e — d-|-1). 

Examples 13.2113.31 and 13.41 suggest that (j3.13() is an equality. 

Conjecture 3.7. The dimension of EC 2 ,d equals d^ — {e — d+ 1). 

We do not know the dimension of 7^2,d = Hfc kernel)^^^]). In our examples, we 
saw that dim(/tr 2 ,d) = 1, 0,3,10 for d = 2, 3,4, 5 respectively. It would be desirable 
to better understand the common kernel Kn,d for arbitrary n and d: 

Problem 3.8. Find the dimension of the space Kn,d in Id. 771) . 


Another problem is to understand the diagonal of EC„^d in the embedding (13.31) . 
This diagonal parametrizes simultaneous eigenconfigurations, arising from special 
tensors A whose d maps ..., all have the same fixed point locus in P”“^. 
Symmetric tensors A have this property, and the issue is to characterize all others. 

Example 3.9 (n = 2). The diagonal of EC 2 ,d is computed by setting = 
= • • • = in the prime ideal described in Theorem l3.61 If d is odd, then there 
is no constraint, by Theorem 12.71 However, for d even, the diagonal of EC 2 ,d is 
interesting. For instance, for d = 2, equating the rows in (13.21) gives two components 


2 X 2-minors of 


do di d2 
02 ~di do 


= (do + d2) ("I (do - d2, ai). 


The hrst component is the known case of symmetric 2 x 2-matrices. The second 
component is a point in Sym 2 (P^), namely the binary form s'^ -\-t^ = {s — it){s + it). 

This is the simultaneous eigenconfiguration of any matrix A = i ^ ^ j with 6^0. 


4. The eigendiscriminant 


The d-dimensional tensors of format nxnx ■ ■ ■ xn represent points in a projec¬ 
tive space P'^ where TV = — 1. For a generic tensor A G P^, the 7*^ eigenconfig¬ 

uration, in the sense of (12.21) . consists of a finite set of reduced points in P”“^. We 
know from Theorem that the number of these points equals 


n—1 

p{n,d,i) = ^(d-l)Y 
2=0 


ifin—l—i 


In this section we study the set An^d,e of all tensors A for which the eigencon¬ 
figuration consists of fewer than p(n, d, £) points or is not zero-dimensional. This 
set is a subvariety of P^, called the 7**^ eigendiscriminant. We also abbreviate 


1 


(4.1) 7(n,d,£) = > (-1) 


7=2 


\n-l-l-J 


O'-i) 


3 


.k—0 


(-1)^ 


n 

J-k, 


j(d—1) — ki — 1 
n — 1 


The following is our main result in this section; 
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Theorem 4.1. The eigendiscriminant is an irreducible hypersurface with 
(4.2) degree(A„,d_^) = 2-f{n,d,£) + 2p{n,d,£) - 2. 

We identify A.n^d,i with the unique (up to sign) irreducible polynomial with 
integer coefficients in the n‘^ unknowns aqia-'-id that vanishes on this hypersurface. 
From now on we use the term eigendisciminant to refer to the polynomial An^d,e- 

The case of most interest is £ = 1, which pertains to the eigenconfiguration of 
a tensor in the usual sense of dOl). For that case, we write = An^d,i for the 
eigendiscriminant, and the formula for its degree can be simplified as follows: 

Corollary 4.2. The eigendiscriminant is a homogeneous polynomial of degree 
degree(A„,d) = n(n - l)(d - 1)”“\ 

The following proof is due to Manuel Kauers. We are grateful for his help. 


Proof. We set £ = 1 in the expression (14.21) . Our claim is equivalent to 


(4.3) 


7 (n,d, 1) = 


(d-1) 


(d-l)--l 


d-2 


1 . 


We abbreviate the innermost summand in (14.11) as 


Sn.dj(fc) = (-1)'‘ 


n — 1 


^ - \j-k, 

Using Gosper’s algorithm [PWZl Chapter 5], we find the multiple 

in.djffcj nj{d-2) 

It can now be checked by hand that this satisfies 

^n^dpik “t“ 1) Sji^dpik^ — 'Sn,dj'(^)- 

Summing over the range fc = 0,..., j — 1 and simplifying expressions lead to 


(4.4) = ^n,d,i(j + l)-5'n,cij(0) = 

This is valid for all j > 1. 

Next we introduce the expression 


d—1/n — 1\/ j(d — 1) — 1 


d — 2 \ j J \ n — 1 


(4.5) 


A(n,d) = ^(-1)^ 
1=0 


n'\ fdj — 1 


(J-I). 


Consider (j — 1) as a polynomial in j of degree n + 1. In the binomial basis. 


dj - 1 
n 


U - 1) = (r + l)d' 


J 

n + ly 

+ lower degree terms. 


+ (d - l)d’ 


n-l I ri+ 1 


-d' 


Recall from [GKPl page 190] that 


1=0 


n] fj 


jj \k. 


{-ly 

0 


if fc = n, 
otherwise. 
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-d^ 


This implies 

A{n,d) = l)d" 

Combining this identity with (14.1|) . (14.41) and (14.51) . we now derive 


2 


'y{n,d, 1 ) = (- 1 )" 


= (- 1 )’ 


d-l 


d-2 

-id-1 

d-2 


A{n — 1, d — 1) + 


-1 
n — 1 


(-i)"-M(d-2)(d-ir 


d-2 
n 
2 


{d-2){d-l) 


n—2 


n 

2 

n—1 




d-l 


d -2 


(d-ir- 


(d- 1 )” - (d- 1 ) 
d -2 ' 



This equals the desired expression for 7 ( 71 , d, 1) on the right hand side of (14.31) . □ 


The proof of Theorem 14.11 involves some algebraic geometry and will be pre¬ 
sented later in this section. We first discuss a few examples to illustrate 

Example 4.3 (d = 2). The eigendiscriminant of an n x n-matrix A = {uij) is 
the discriminant of its characteristic polynomial. In symbols, 

= discrA(det( A — A • Id„)). 

This is a homogeneous polynomial of degree n(n — 1) in the matrix entries Uy. For 
instance, for a 3 x 3-matrix, the eigendiscriminant is a polynomial with 144 terms: 

A3,2 = 0,110-22 ~ 2 aiia 22 U 33 + 023032 + O11O33 ~ 2a^]^ai2a2ia22 -f • • • -f 023^32^33- 

This polynomial vanishes whenever two of the eigenvalues of A coincide. 

There is a beautiful theory behind A „,2 in the case when A is real symmetric, so 
the eigenconfiguration is defined over R. The resulting symmetric eigendiscriminant 
is a nonnegative polynomial of degree n(n — 1) in the matrix entries. Its real 

variety has codimension 2 and degree , and its determinantal representation 
governs expressions of A „,2 as a sum of squares of polynomials of degree ( 2 ). For 
further reading on this topic see [Stul Section 7.5] and the references given there. 

Example 4.4 (n = 2). The eigendiscriminant of a d-dimensional tensor of 
format 2 x 2 x • • • x 2 is the discriminant of the associated binary form in ( 12 .7L i.e. 

A2,d = disc(:j,y)(?/• 7/11(0;, y) - a; •'(/'2(a;,y))- 

This is a homogeneous polynomial of degree 2d — 2 in the 2'^ tensor entries 

Example 4.5 (n = d = 3). The eigendiscriminant A 3,3 of a 3x3x3-tensor 
A = {oijk) is a homogeneous polynomial of degree 24 in the 27 entries Oijk- If we 
specialize ^ to a symmetric tensor, corresponding to a ternary cubic 

<p{x,y,z) = csoox^ + C 2 iox‘^y + C 2 oix'^z-\ -hcoos-z^, 

then A 3,3 remains irreducible. The resulting irreducible polynomial of degree 24 in 
the ten coefficients Cyfe is the eigendiscriminant of a tenary cubic. At present we 
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do not know an explicit formula for but it is fun to explore specializations of 
the eigendiscriminant. For instance, if (j) = {2x + y){2x + z){2y + z)+u ■ xyz then 

A 3,3 = 1611^'* + 2304u^^ + 15278411^^ + 6097536ii^^ + 159761808u^'’ + 2779161840ii^® + 29727588168u^® 

+ 124641852624u'^’' - 123407858901611^® - 18314627517360u'^® - 8929524942432ii^<‘ + 120093304792564811^® 
+3722203539791685U®® - 63418425922462464u®® - 257381788882972176ii®® + 2676970903961440800u® 

+ 792765511483628649611® - 89013482239908955392ti'^ - 13934355026171012352u® 

+ 172925035637155679232011® - 515922232490119293004811“ - 11838757458480721920u® 
+2825545664173411698278411® - 56809371779894977339392ii + 37304830510913780269056, 

and if (j) = u ■ x^ + v ■ y^ + w ■ z^ + xyz then Aa,a is the square of polynomial 

531441ii“tj“iii“ - 70858811®tJ®iD® - 708588 ii®ii®id® - 708588ii® ii®id® + 1062882ii“ii“iii® + 1062882ii“tJ®iii“ 

+ 106288211®u“iii“ - 1810836ti®u®iD® - 177147ti“u“ + 39366ii“D®m® + 3936611®!)“™® - 177147 ii“id“ 
+3936611®!)®™“ - 1771471)“™“ + 314928ii®d®™ + 314928ii®d™® + 314928iid®™®-244944ii®d®™®-46656ii®d™ 
-46656111)® ™-46656iid™®+2332811®!)®+2332811® ™®+23328d®™®+6912iid™- 230411® - 2304t)® - 2304™® + 256. 

We now embark towards the proof of Theorem 14.11 Let i4 and B be generic 
tensors of the same format, and write (i/'i ,... ,'ijjn) and (wi ,... ,ujn) for the vectors 
of degree d—1 polynomials that represent the corresponding maps 
Let C denote the subvariety of P”“i defined by the determinantal constraints 

M(x) l/>2(x) ••• V'nWN 

(4.6) rank a;i(x) u; 2 (x) ••• uJn(x) < 2. 

\ T-f- T-f- ... ] 

\ ■’'1 -^2 •'-ra / 

Since the ipi and ujj are generic, this defines a variety of codimension n — 2. We 
find that C is a curve that is smooth and irreducible, by an argument similar to 
that in the proof of Theorem 12.21 The following lemma is the key to Theorem 14.11 

Lemma 4.6. The expression in o is the genus of the curve C. In symbols, 

genus(C) = j{n,d,£). 

Using this lemma, we now derive the degree of the eigendiscriminant. 

Proof of Theorem 14.11 We define a map /r : C —?> P^ as follows. For any 
point X on the curve C, the matrix in ( 14 . 61 ) has rank 2 , so, up to scaling, there exists 
a unique row vector (a, b, c) G that spans the left kernel of that 3 x n-matrix. We 
define the image of x G C to be the point /r(x) = (a : 6) on the projective line P^. 
This condition means that x is an eigenpoint of the tensor a A + bB. Conversely, 
for any (a : 6) G P^, the fiber : b) consists precisely of the eigenpoints of 

aA + bB. Hence, since A and B are generic, the generic fiber is finite and reduced 
of cardinality p{n, d, £). In other words, /r : C —)• P^ is a map of degree p{n, d, £). 

We restrict the eigendiscriminant to our P^ of tensors. The resulting binary 
form /S.n^d,i{aA + hB) is squarefree, and its degree is the left hand side in (14.21) . The 
points {a : b) G P^ where An^d.iiaA + bB) = 0 are the branch points of the map p. 
The corresponding multiplicity-two eigenpoints x form the ramification divisor on 
C. The number of branch points of p is the degree of the eigendiscriminants An^d,e- 
The Riemann-Hurwitz Formula IM Exercise 8.36] states that the number of 
branch points of the map /i : C —>• P^ is 2 ■ genus(C) -I- 2 • degree(/r) — 2. By the first 
paragraph, and by Lemma [4.61 this expression is the right hand side of (14.21) . □ 

Our proof of Lemma 14.61 is fairly complicated, and we decided not to include 
it here. It is based on resolutions of vector bundles, like those seen in the proof of 
Theorem 12.21 We plan to develop this further and publish it in a later paper on 
discriminants arising from maximal minors of matrices with more than two rows. 

What we shall do instead is to prove an alternative combinatorial formula for 
the genus of C that is equivalent to (14.11) . This does not prove Lemma l476l because we 
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presently do not know a direct argument to show that they are equal. Nevertheless, 
the following discussion is an illustration of useful commutative algebra techniques. 

Instead of Lemma 14.61 The Hilbert polynomial Hc{t) of the curve C equals 

He{t) = degree(C) • t + (1 — genus(C)). 

Recall that C is a linear section of the variety defined by the maximal minors of a 
3 X n-matrix whose rows are homogeneous of degrees d — 1, d — 1 and That variety 
is Cohen-Macaulay. We shall compute the Hilbert polynomial of the coordinate ring 
of C from its graded minimal free resolution over S = K[xi ,..., Xn]- 

Consider the S'-linear map from F = S'®" to G = S(d — 1)®^ © S(£) given by 

V'2(x) ••• l/’„(x)^ 

(4.7) a= I wi(x) W 2 (x) ••• w„(x) 


By [Eli Section A2H], the corresponding Eagon-Northcott complex EN(Qf) equals 


0 ^ Sym„_3(G^)©/\F ^ Sym„_4(G'^)© /\F 


G''^/\F /\F^ /\G, 


where Symj(G^) is the **** symmetric power of G^ and /\* F is the z**' exterior power 
of F. We compute the Hilbert polynomial HM{t) of each module M in EN(a). 

Since C has codimension n — 2 = rankF — rankG+1, the complex EN(a) © 
S{—2d — ^ + 2) is a free resolution of the coordinate ring of C. In particular, 

n 

(4.8) Hc{t) = H/^3G(<-2d-£ + 2) + ^(-l)^id£;.(t-2d-£ + 2), 

i=3 

where Ej = Sym■_ 3 (G^) © /\^ F. Since G^ = S{—d + 1)®^ i 


) S{—£) and F = 


3-3 


Sym^._ 3 (G'^) = 0S((j-fc-3)(-d+l)-H)®^-'=+' and /\e = S®("). 




Their tensor product is the term in EN(a). As a graded 5-module, it equals 
i-3 

Ej = 0S((j-/c-3)(-d+l)-H)®^^“''“^^(^). 


fc =0 


The shifted Hilbert series of this module is the summand on the right of (14.81) : 

3-3 


HE.{t-2d-£ + 2) = ^(j-fc-2)( . 

fe=0 


t -h (n—1) -h {j—k—l){l—d) — 
n — 1 


■ 


We conclude that the Hilbert polynomial Hc{t) of the curve C equals 

n 3-3 

3=3 


t+n—1 
n — 1 


k=0 


t + (n-1) + {j-k-l){l-d) - (fc+l)A 
n — 1 / 


The genus of C is obtained by substituting t = 0 and subtracting the result from 1: 

n j-3 

7(n,d,£) = ^(-l)^“^^(j-fc-2)( 

j—3 k—0 


(n-l) + (j-fc-l)(l-d) - {k+l)£\ 
n — I 


■ 


This formula is equivalent to grj. 


□ 
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Example 4.7 (n = 4). The last formula seen above specializes to 


7 ( 4 , (i,£) = 4 


3 + 2(l-d) 
3 


3 + 3(l-d) 
3 


3 + 2(l-d) - 2£ 
3 


while the genus formula in (14.11) states 


7 ( 4 , d,£) = 


+2 


4(3(d-3l)-l) _ _ ^3(d-lK3f-l^ 


Both of these evaluate to the cubic polynomial 

7(4, d, 1) = + 3d£^ + t - 21d^ - Ud£ - + 27d + 9£ - 10. 

Therefore, by Theorem l4.ll the degree of the £*** eigendiscriminant for n = 4 equals 

degree(A4.d7) = Sd^ + + 2df + £^ - 12d^ - Ml - + I5d +M-6. 

For 1 = 1 , this factorizes as promised in Corollary 14.21 degree(A 4 ^d 4 ) = 12(ci— 1)^. 


5. Seven points in the plane 

Our study had been motivated by the desire to find a geometric characterization 
of eigenconfigurations among all finite subsets of P”“^. The solution for n = 2 was 
presented in Theorem l2.7l However, the relevant geometry is more difficult in higher 
dimensions. In this section we take some steps towards a characterization for n = 3. 
The eigencoiffiguration of a general tensor A in consists of d^—d+1 points in 

P^. So, our question can be phrased like this: given a configuration Z € (P^)'^ -d+i^ 
decide whether it is an eigenconfiguration. If yes, construct a corresponding tensor 
A € and decide whether A can be chosen to be symmetric. 

The first interesting case is d = n = 3. Here the following result holds. 

Theorem 5.1. A configuration of seven points in P^ is the eigenconfiguration 
of a 3 X a X S-tensor if and only if no six of the seven points lie on a conic. 

The only-if part of this theorem appears also in mi Proposition 2.1], where 
Ottaviani and Sernesi studied the degree 54 hypersurface of all Luroth quartics 
in P^. We note that part (i) in m Proposition 2.1] is not quite correct. A 
counterexample is the configuration Z consisting of four points on a line and three 
other general points. It is precisely this gap that makes our proof of Theorem 15.11 
a bit lengthy. 

This proof will be presented later in this section. Example 11.11 shows that some 
triples among the seven eigenpoints in P^ can be collinear. Another interesting 
point is that being an eigenconfiguration is not a closed condition. For a general d 
it makes sense to pass to the Zariski closure. We define Eig^ to be the closure in 
(p2)‘^ -d+i qj: eigenconfigurations. Readers from algebraic geometry 

may prefer unlabeled configurations, and they would take the closure in the Chow 
variety Sym^ 2 _(i_|_ 4 (P^) or in the Hilbert scheme Hilb£; 2 _(j+ 4 (P^). For simplicity of 
exposition, we work in the space of labeled point configurations. We also consider 
the variety of symmetric eigenconfigurations, denoted Eig^^ g^j^,,. This is the Zariski 
closure in (P^)'^ -d+i gg^ eigenconfigurations of ternary forms f of degree d. 

Towards the end of this section we examine the dimensions of Eig^ and Eig^ 
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We begin by approaching our problem with a pinch of commutative algebra. 
Let Z G -d+i write Iz for the ideal of all polynomials in S = K[x,y^z\ 

that vanish at all points in the configuration Z. This homogeneous radical ideal is 
Cohen-Macaulay because it has a free resolution of length 1 (see, for example, [Eisl 
Proposition 3.1]). By the Hilbert-Burch Theorem, the minimal free resolution of 
Iz has the form 

0 ^ A S®'^ ^ Iz ^ 0. 

The m x (m—l)-matrix <I> is the Hilbert-Burch matrix of Z. The minimal free reso¬ 
lution of Iz is unique up to change of bases in the graded S'-modules. In that sense, 
we write <Pz The ideal Iz is generated by the maximal minors of <I>z- The 

following proposition is due to Ottaviani and Sernesi (see mi Proposition 2.1]). 

Proposition 5.2. Let Z be a configuration in (p2)rf^-d+i^ Then Z is the 
eigenconfiguration of a tensor if and only if its Hilbert-Burch matrix has the form 

/ Li Fi\ 

(5.1) = \ L2 F2 \, 

\ L3 F3 ) 

where Li,L2,L3 are linear forms that are linearly independent over K. 

This statement makes sense because the condition on <I>z is invariant under row 
operations over K. The ternary forms ^ 1 ,^ 2 , P’s must all have the same degree, 
and the hypothesis on Z ensures that this common degree is d — 1. 

Proof. We start with the only-if direction. Suppose that Z is an eigenconfig- 
uration. Then there exist ternary forms 4’i,4’2,4>3 of degree d — 1 such that Z is 
defined set-theoretically by the 2 x 2-minors of 

(5 2) ( y ^ ] 

yPiix,y,z) '4j2{x,y,z) tp3{x,y,z)J 

The ideal generated by these minors is Cohen-Macaulay of codimension 2 and its 
degree equals the cardinality of Z. This implies that this ideal coincides with Iz- 
The Hilbert-Burch Theorem ensures that the transpose of (|5.2I1 equals I^z- Since 
x,y,z are linearly independent, we see that <I>z has the form required in (|5.1ll . 

For the converse, suppose that the Hilbert-Burch matrix of Z has size 3 x 2 as 
in EID with Li,L 2 ,L 3 linearly independent. By performing row operations over 
K, we can replace Li, T 2 , L3 by x, y, z. This means that the transpose of <I>z is (15.21) 
for some V'l j '02, V '3 • Any such triple of ternary forms of degree d — 1 arises from 
some tensor A G By construction, Z is the eigenconfiguration of A. □ 

Proposition 15.21 translates into an algorithm for testing whether a given Z G 
(p 2 )d -d-i-i jg eigenconfiguration. The algorithm starts by computing the ideal 

Iz = Pi {x[ 3 - ya, xj - za, yj - z/3). 

{a:P:-y)GZ 

This ideal must have three minimal generators of degree d; otherwise Z is not an 
eigenconfiguration. If Iz has three generators, then we compute the two syzygies. 
They must have degrees 1 and d — 1, so the minimal free resolution of Iz looks like 

0 ^ S{-d - 1) © S'(-2d -t 1) A S'(-d)®3 ^ Iz ^ 0- 
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At this point we examine the matrix <1>. If the linear entries Li,L 2 ,L 3 in the left 
column are linearly dependent, then Z is not an eigenconfiguration. Otherwise we 
perform row operations so that <P^ looks like (15.21) . The last step is to pick a tensor 
A € that gives rise to the ternary forms , ■i/' 2 , V's in the second row of . 

The remaining task is to find a geometric interpretation of the criterion in 
ProDOsition l5.2l This was given for d = 3 in the result whose proof we now present. 

Proof of Theorem 15.11 Fix a configuration Z G (P^)^. Our claim states 
that the Hilbert-Burch matrix (!>z has format 3 x 2 as in with Li,L 2 , L 3 

linearly independent, if and only if no six of the points in Z lie on a conic. 

We begin with the only-if direction. Take p G Z such that Z\{p} lies on a 
conic C in P^. Fix linear forms Li and L 2 that cut out p. The cubics CLi and 
CL 2 vanish on Z. By Proposition 15.21 we have Iz = {CLi,CL 2 ,F) where F is 
another cubic. Since F vanishes at p, there exist quadrics Qi and Q 2 such that 
F = Q 2 L 1 — QiL 2 . The generators of the ideal Iz are the 2 x 2-minors of 

/ Li Qi \ 

F = \ L 2 <52 I ■ 

\ 0 c J 

This means that is a Bilbert-Burch matrix <Pz for Z. Bowever, by Proposition l5.21 
the left column in any Fz must consist of linearly independent linear forms. This 
is a contradiction, which completes the proof of the only-if direction. 

We now establish the if direction. Fix any configuration Z G (P^)^ of seven 
points that do not lie on a conic. We first prove that the minimal free resolution 
of Iz has the following form, where c is either 0 or 1 : 

(5.3) 0 ^ © S{-5) -G 5'(-3)®3 © ^(-d)®^ ^ 0,. 

By the Bilbert-Burch Theorem, the resolution oi Iz equals 

t t+i 

where t, oi,..., Oj+i, bi,... ,bt G N with oi > • • • > at+i and bi > ■ ■ ■ > bt- We 
abbreviate et =bi — at and fi = bi — a^+i for z G {1,..., t}. These invariants satisfy 

(i) fi ^ Bi, Ci+i, 

(ii) Bi,fi > 1 . 

Furthermore, Eisenbud shows in |Eisl Proposition 3.8] that 
(hi) © = +x;*=i/j- 

There exist 3 linearly independent cubics that vanish on the seven points in Z. By 
[Eisl Corollary 3.9], the ideal Iz has either 3 or 4 minimal generators, so t G {2, 3}. 
Suppose t = 2. Then oi = 02 = 03 = 3, and it follows from (hi) that 

/i + /2 = 3 and ei + /2 = 3. 

So, by (i) and (ii), we obtain ei = /i = 1 and /2 = 2. This implies 61 = 5 and 
62 = 4. Therefore, Iz has a minimal free resolution of type (15.31) with c = 0. 

Next, suppose t = 3. Then oi > 02 = 03 = 04 = 3. From (hi) we now get 

/i + /2 + /a = «! 

ei + /2 + /3 = 3 

ei + €2 + h = 3. 
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By (ii), ei = 62 = /2 = /s = 1- Therefore, 62 = 02 + 62 = 4 and 63 = 02 + /2 = 4. 
Corollary 3.10 in |Eisj says that 

eifj = 61 (/i + /2 + /a) + 62(/2 + /a) + 63/3 = oi + 3 = degZ = 7. 

i<j 

Hence oi = 4, 61 = 5, and Iz has a minimal free resolution of type (15.31) with c = 1. 

To complete the proof, we now assume that no six points of Z lie on a conic. 
In particular, no conic contains Z, so the minimal free resolution oi Iz equals (15.3L 
with 6 € {0,1}. Suppose that c = 1. The Hilbert-Burch matrix must be 




Too 

Toi 

Qo 

Tio 

Til 

Qi 

T20 

T21 

Q2 

0 

0 

L 


with L, Lij are linear and Qk are quadrics. 


Then Iz = {LQq, LQ[, LQ 2 , Q), where 


Qo = 


Too 

Tio 


Toi 

Til 


,Q'i = 


Too 

T 20 


Toi 

T 21 


Q2 


Tio 

T 20 


Til 

T 21 ’ 


Q 


Too Toi Qo 
Tio Lii Qi 
T20 T21 Q2 


The ideal generated by L and Q contains Iz- The intersection of the curves 
{L = 0} and {Q = 0} is contained in Z. Note that these curves share no positive¬ 
dimensional component, since Z is zero-dimensional. Thus {L = Q = 0} consists of 
four points. Let L' be a linear form vanishing on two of the three other points. Then 
the conic {LL' = 0} contains six points of Z, which contradicts our assumption. 
Hence, c = 0. The resolution (15.31) tells us that the Hilbert-Burch matrix equals 

/ Ti Qi \ 

= I T 2 Q 2 I, 

\ T3 Qs J 

with linear forms Li and conics Qj. If Li,L 2 ,T 3 were linearly dependent then we 
can take L 3 = 0. So, the conic Q 3 contains the six points in Z\{Li = L 2 = 0}. 
Consequently, the linear forms Li,L 2 ,L^ must be linearly independent. Proposi¬ 
tion O now implies that Z is the eigenconliguration of some 3 x 3 x 3-tensor. □ 


After taking the Zariski closure, we have Eigg = (P^)^. We shall now discuss the 
subvariety Eigg of those eigenconfigurations that come from symmetric tensors. 
Consider the three quadrics in the second row of (15.21) . We write these as 

ipiipL) = aix^ + a2xy + asxz + + a^yz + , 

V'2(x) = bix'^ + b2xy + b^xz + b^y"^ + b^yz + b^z'^, 

V'sW = cix"^ + C2xy + c^xz + c^y"^ + c^yz + cez"^. 

We shall characterize the case of symmetric tensors in terms of these coefficients. 

Proposition 5.3. The variety Eigg ^yjj^ is irreducible of dimension 9 in (P^)^. 
An eigenconfiguration Z comes from a symmetric tensor (j) as in if and only if 

0-5 ~ b 3 = 1>3 — 62 = 204 ~ 206 ~ ^2 + 63 = 265 ~ 26i — 65 -|-02 = 2ci — 2c4 -03-1-65 = 0. 


Proof. There exists a symmetric tensor with eigenconliguration Z if and only 
if there exist a cubic (j) and a linear form L = uix + U 2 y -f u^z such that 


(5.4) 


IT , T , T ^4’ 

tpi+ Lx = —, Ip2 + Ly = —, ili 3 + Lz=—. 

ox ay oz 
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We eliminate the cubic (j) from this system by taking crosswise partial derivatives: 
dipi dL dilj 2 dL 9'(/>2 dL dtps dL 
^ ^ ^ ’ ■ ■ ■ ’ ^ ^ ^ ^ ^ ■ 

This is a system of linear equations in the 21 unknowns oi,..., oe, &i,..., 651 ci, ..., 
Cq,Ui,U2,us. By eliminating the last three unknowns Ui,U2,U3 from that system, 
we arrive at the hve linearly independent equations in a^, bj, Ck stated above. □ 

Proposition 15.31 translates into an algorithm for testing whether a given config¬ 
uration Z is the eigenconfiguration of a ternary cubic (p. Namely, we compute the 
syzygies of Iz, we check that the Hilbert-Burch matrix has the form (15.21) . and then 
we check the five linear equations. If these hold then (p is found by solving dSH). 

While the equations in Proposition l5.3l are linear, we did not succeed in comput¬ 
ing the prime ideal of Eigg in the homogeneous coordinate ring of (P^)^. This 
is a challenging elimination problem. Some insight can be gained by intersecting 
Eigg gym with natural subfamilies of (P^)^. For instance, assume that Z contains 
the three coordinate points, so we restrict to the subspace (P^)^ defined by 

Z = {(1:0:0), (0:1:0), (0:0:1), (tti: 02 : 03 ), {Pi-h'Pti), ( 71 : 72 : 73 ), (^i:^2:^3)}- 

At this point it is important to recall that our problem is not projectively invariant. 

Theorem 5.4. The variety Eigg gym H (P^)'^ is three-dimensional, and it repre¬ 
sents the eigenconfigurations Z of the ternary cubics in the Hesse family 

(5.5) (p = ax^ -\- by^ -\- cz^ -\- Mxyz. 

If a, b, c, d are real then the eigenconfiguration Z contains at least five real points. 

Proof. A ternary cubic p = J2i+j+k=3 has (1 : 0 : 0) as an eigen- 

point of p if and only if C 210 = C 201 = 0. Likewise, (0 : 1 : 0) is an eigenpoint if and 
only if C 120 = C 021 = 0, and (0 : 0 : 1) is an eigenpoint if and only if C 102 = C 012 = 0. 
Hence the eigenconfiguration of p contains all three coordinate points if and only 
if p is in the Hesse family (15.51) . Since Eigg gym has codimension 5 in (P^)^, the 
intersection Eigg gym H (P^)^ has codimension < 5, so its dimension is > 3. The 
Hesse family is 3-dimensional, and so we conclude that dim(Eigg_gym n (P^)"^) = 3. 

The four other eigenpoints of (15.51) are {dxix^ — 1) : x(«X “ c) : d(x^ — 1)), 
where x runs over the zeros of the polynomial 

(5.6) d{a^ — d^)x^ — a{ab-\-cd)^ -\- 2{abc-\-d^)'x^ — c{bc-\-ad)x + d{(? — d?). 

We claim that this quartic polynomial has at least two real roots for all a, b,c,d € M. 

Inside the projective space P^ of quartics /(x) = fc 4 X^ + ^ 3 X^ + ^ 2 X^ + fciX + ^o, 
the family (15.61) is contained in the hypersurface defined by the quadric 

tS = 2/C4/C2 -)- k 2 — d/Cg/Cg -)- -\- 2/Cq/C2. 

The discriminant of /(x) defines a hypersurface of degree 6 in P"^. One of the open 
regions in the complement of the discriminant consists of quartics /(x) with no 
real roots. In polynomial optimization (cf. |BPTl Lemma 3.3]) one represents this 
region by a formula of the following form, where k is a new indeterminate: 

/fco fci/2 K \ /x^\ 

fix) = (x^ X 1) • ki k2 - 2 k ksl 2 ■ X , 
y K ksl2 ki j \l j 
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The symmetric 3 x 3-matrix is required to be positive definite for some k € K. The 
condition of being positive definite is expressed by the leading principal minors: 

/ fco A:i/2 K \ 

, TZ = det I ki/2 k 2 —2k ^ 3 / 21 . 

\ K kzl2 ki ) 

It remains to be proved that there is no solution (feg, fcr, ^ 2 ) ^ 3 , ^ 4 , S K® to 

V > Q, Q>0, 'R.> 0 and 5 = 0. 

We showed this by computing a sum-of-squares proof, in the sense of |BPTl Chapter 
3]. More precisely, using the software SDSTools [SOSj . we found explicit polyno¬ 
mials p, q,r, s G M[/co, ki, ^ 2 , /ca, ^ 4 , k] with floating point coefficients such that 

p, q, r are sums of squares and pV + qQ + rTZ -|- s5 = —1. 

We are grateful to Cynthia Vinzant for helping us with this computation. □ 

We close this section by returning to tensors in for general d> 3. 


'P = ko, Q = det 


ko ki/2 
ki/2 k2 — 2 k 


Theorem 5.5. Consider the spaees of eigenconfigurations of ternary tensors, 
C Eig, C for d > 3. 

These projective varieties are irreducible, and their dimensions are 

dim(Eig^ gyjjj) = -(d^ -|- 3d) and dim(Eig^) = d^ -I- 2d — 1. 

Proof. First we show dim(Eig^) = d^ -|- 2d — 1. Let W be the set of 2 x 3 
matrices (15.21) . This is a 3 ('^ 2 ^)-dimensional vector space over K. The group 


G = 


1 0 
/ a 


a G isr\{0} and / is a ternary form of degree d — 2 


acts on W by left multiplication. Consider (p,[j G W. It is immediate to see that if 
uj = g ■ ip for some g G G, then the variety defined by the 2 x 2-minors of p equals 
the variety defined by the 2 x 2-minors of uj. The converse also holds because of 
the uniqueness of the Hilbert-Burch matrix. The set W° of elements in W whose 
2 X 3-minors define d^ — d -I- 1 distinct points is an open subset of W. Therefore, 


dim(Eigrf) = dimW7G = dimW°-dimG 


= 3 


d -I-1 
2 


= d^ -b2d- I. 


Next we prove dim(Eig^ 

I / T* -7/ 

U := 


) = ^(d^ + 3d). We introduce the linear subspace 


X y z 

dcfjdx dcj/dy dtfjdz 

The action of the group G on VE does not restrict to U. In fact, we notice that 

I 0 \ / X y z 

f a y d(j)jdx dcf/dy dcf/dz 

if and only if / = 0. Let U° = U Ci W° and consider the subgroup 

I 0 


(p ternary form of degree d > C VE 


G U 


H = 


0 a 


a G K\{0} 


C G. 
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This yields dim(Eig^ gy^) = d\TaU°/H = — 1 = +3d), as desired. Our 

configuration spaces Eig^ and Eig^^ ^yj^j are irreducible varieties because they contain 
the irreducible varieties W° jG and U°/H respectively as dense open subsets. □ 

6. Real eigenvectors and dynamics 

In this section we focus on the real eigenpoints of a tensor A in If A is 

generic then the number of eigenpoints in equals ((d — I)" — l)/(d — 2). Our 
hope is to show that all of them lie in for suitably chosen symmetric tensors 
(j). A second question is how many of these real eigenpoints are robust^ in the sense 
that they are attracting fixed points of the dynamical system Vcj) : Pr~^ —->■ Pr~^. 
Our results will inform future numerical work along the lines of |CDN1 Table 4.12]. 

We begin with a combinatorial construction for the planar case (n = 3). Con¬ 
sider an arrangement A of d distinct lines in PL and let </> be the product of d linear 
forms in x, y, z that define the lines in A. We assume that A is generic in the sense 
that no three lines meet in a point. Equivalently, the matroid of A is a uniform 
rank 3 matroid on d elements. Such an arrangement A has ( 2 ) vertices in P^, and 
these are the singular points of the reducible curve {(j) = 0}. The complement of A 
in Pg has ( 2 ) + 1 connected components, called the regions of A. 

We are interested in the eigenconfiguration of A, by which we mean the eigen- 
configuration of the symmetric tensor ij). Theorem 12.11 gives the expected number 

(6.1) d^—d-l-1 = 1 -I- (d—1) -I- (d—1)^ = 2^^^-1-1 =# vertices-|-# regions. 

The following result shows that this is not just a numerical coincidence. 

Theorem 6.1. A generic arrangement A of d lines otPr has d^ — d+1 complex 
eigenpoints and they are all real. In addition to the ( 2 ) vertices, which are singular 
eigenpoints, each of the ( 2 ) + 1 regions of A contains precisely one real eigenpoint. 

Proof. The singular locus of the curve {(j) = 0} consists of the vertices of the 
arrangement A. These are the eigenpoints with eigenvalue 0. Their number is ( 2 ). 

Let Li, L 2 , ■. ■, Ldhe the linear forms that define the lines, so </) = L 1 L 2 • ■ ■ Ld- 
Consider the following optimization problem on the unit 2-sphere: 

d 

Maximize log I(/)(x)I = ^ log |Li(x, y, z)] subject to -I- -f = 1. 

The objective function takes the value —00 on the d great circles corresponding to 
A. On each region of A, the objective function takes values in R, and is strictly 
concave. Hence there exists a unique local maximum u* = (x*, y*, z*) in the interior 
of each region. Such a maximum u* is a critical point of the restriction of (/'(x) to 
the unit 2-sphere. The Lagrange multiplier conditions state that the vector u* is 
parallel to the gradient of at u*. This means that u* is an eigenvector of (j), and 
hence the pair ±u* defines a real eigenpoint of (j) in the given region of Pg. 

We proved that each of the ( 2 ) + 1 regions of A contains one eigenpoint. In 
addition, we have the ( 2 ) vertices. By Theorem 12.11 the total number of isolated 
complex eigenpoints cannot exceed 2 ( 2 ) -f 1. This means that there are no eigen¬ 
points in Pp other than those already found. This completes the proof. □ 
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We note that the line arrangement can be perturbed to a situation where the 
map S/cj) ■ regular, i.e. none of the eigenvectors has eigenvalue zero. 

Corollary 6.2. There exists a smooth curve of degree d in the real projeetive 
plane whose complex eigenconfiguration consists of — d + 1 real points. 

Proof. The eigenconfiguration of (/) = L 1 L 2 • ■ ■ Lj, is 0-dimensional, reduced, 
and defined over K. By the Implicit Function Theorem, these properties are pre¬ 
served when (j) gets perturbed to a generic ternary form that is close to (f. □ 


It is interesting to see what happens when the matroid of A is not uniform. Here 
the eigenconfiguration is not reduced. It arises from Theorem 16. 1 1 bv degeneration. 

Example 6.3. Let d = 6 and take A to be the line arrangement defined by 
cj) = X ■ y z ■ {x - y) ■ {x - z) ■ {y - z). 

This is the reflection arrangement of type A 4 . Its eigenscheme is non-reduced. 
Each of the 12 regions contains one eigenpoint as before, and the simple vertices 
(1 : 1 : 0), (1 : 0 : 1), and (0:1:1) are eigenpoints of multiplicity one. However, 
each of the triple points (1 : 0 : 0), (0 : 1 : 0), (0 : 0 : 1), (1 : 1 : I) is an eigenpoint of 
multiplicity 4. This makes sense geometrically: in a nearby generic arrangement, 
such a vertex splits into three vertices and one new region. We note that the scheme 
structure at the eigenpoint (1 : 0 : 0) is given by the primary ideal { 2 yz — z'^, y^ — z^). 

The concavity argument concerning the optimization problem in the proof of 
Theorem 16.II works in arbitrary dimensions, and we record this as a corollary. 


Corollary 6.4. Each of the open regions of an arrangement A of d hyper¬ 
planes in Pg~^ contains precisely one real eigenpoint of A. The number of regions is 


( 6 . 2 ) 


n —1 

E 

2=0 


d — 1 


Proof. The first part has the same proof as the one for n = 3 given above. 
The formula for the number of regions can be found in [stil Proposition 2.4]. □ 


Theorem IQ is restricted to n = 3 because hyperplane arrangements are sin¬ 
gular in codimension 1. Hence the eigenconfiguration of a product of linear forms 
in n > 4 variables has components of dimension n — 3 in We conjecture that 

a fully real eigenconfiguration can be constructed in the vicinity of such a tensor. 


Conjecture 6.5. Let (f be any product of d nonzero linear forms in Sym]^(R”). 
Every open neighborhood of (j) in Sym^(]R"') contains a symmetric tensor such 
that all ((d — 1)” — l)/(d — 2) complex eigenpoints of are real. 

This optimistic conjecture is illustrated by the following variant ofExample ll.il 

Example 6.6 (n = d = 4). The classical Cremona transformation in is V()> 
where p = xyzw is the product of the coordinates. The eigenconfiguration of 4> 
consists of eight points, one for each sign region in and the six coordinate lines. 
The expected number (12.11) of complex eigenpoints is 40. Consider the perturbation 

(pe = xyzw -f e(5a;^ -I- Ax^y — 2x‘^y‘^ — 8 xy^ ly'^ -I- -I- 2x‘^yz -|- 2xy‘^z 

-\- 2 y^z - 6 x^z^ Qxyz"^ -|- ly'^z^ - 8 xz^ 3yz^ Sz"^ - 8 x^w -|- 2 x'^yw 
—3xy^w by^w 8 x‘^zw — 3y^zw — bxz^w — lOyz^w 8 z^w — bx^w^ 
— 6 xyw^ — 3y^w^ — 6 xzw^ 3yzw^ 3xw^ -\- 3yw^ — 4zw^ -|- 3w^). 
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All 40 complex eigenpoints of this tensor are real, so Conjecture 16.51 holds for (f). 

Remark 6.7. Conjecture 16.51 is true for n < 3. For n = 3 this follows from 
Theorem l6.ll we can take tpe to be any perturbation of the given line arrangement cj). 
For n = 2 we take (j)^ = (j) itself because of the following fact: if a binary form (j){x, y) 
is real-rooted then also ~ yfy is real-rooted. This follows from Corollarv l6.4l 

We put the lid on this paper with a brief discussion of the dynamical system 
i/) : --■» associated with a tensor A € Iterating this map is 

known as the tensor power method, and it is used as a tool in tensor decomposition 
[AG]. This generalizes the power method of numerical linear algebra for computing 
the eigenvectors of a matrix A G (R.")®^. One starts with some unit vector v and 
repeatedly applies the map v For generic inputs A and v, this iteration 

converges to the eigenvector corresponding to the largest absolute eigenvalue. 

Suppose that u G Pr~^ is an eigenpoint of a given tensor A G We say 

that u is a robust eigenpoint if there exists an open neighborhood 7/ of u in P^” 
such that, for all starting vectors v G W, the iteration of the map ijj converges to u. 

Example 6.8 (Odeco Tensors). A symmetric tensor is orthogonally decompos¬ 
able (this was abbreviated to odeco by Robeva |Robj 1 if it has the form 

(j) = 

where ai,..., G M and {vi, ..., v^} is an orthogonal basis of Following 
[AG], the robust eigenpoints of an odeco tensor are the basis vectors v^, and they 
can be computed using the tensor power method. Up to an appropriate change 
of coordinates, the odeco tensors are the Fermat polynomials in (j2.4|) . The robust 
eigenpoints oi (j) = xf x‘}^ are the coordinate points ei,..., e„. The region 

of attraction of the eigenpoint under the iteration of the map is the set 
of all points in whose coordinate is largest in absolute value. 

Odeco tensors for n = d = 3 have three robust eigenvalues. At present we do not 
know any ternary cubic (j) with more than three robust eigenpoints. Theorem 16.11 
might suggest that products of linear forms are good candidates. However, we 
ran experiments with random triples of lines in P^, and we observed that the 
number of robust eigenvalues is usually one and occasionally zero. We never found 
a factorizable ternary cubic cj) with two or more robust eigenpoints. The Cremona 
map in Example 11.11 shows that cj) = xyz is a cubic with zero robust eigenpoints. 
Here is a similar example that points to the connection with frame theory in fORS] . 

Example 6.9 {n = d = 3). We consider the factorizable ternary cubic 

(6.3) (j> = {2x2y — z){2x — y2z){—x2y2z). 

This equals the frame decomposable tensor seen in [ORSl Examples 1.1 and 5.2]: 

(6.4) (j> = F((_5a;_|_y_|_2;)3 + jx-5y-\-z)^ {x-\-y-5z)^ -|- (3a;-I-3y-I-3z)^). 

Its gradient map P^ —-> P^ is given by 

/—4a:^ -I- 4:xy -j- Axz 2y^ -\-yz 2z^\ 

S/cj) = 3 j 2a:^ -I- 'ixy xz — 4y^ -|- Ayz 2z^ j . 

y 2a:^ -I- a;y -I- 4x2 -I- 2y^ -|- Ayz — 4z^ j 
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This has four fixed points and three singular points, for a total of seven eigenpoints: 

(l:l:-5), (1:-5:1), (-5:1:1), (3:3:3), 

( 2 : 2 : - 1 ), (2 :-1 : 2 ), (- 1 : 2 : 2 ). 

Note that the pairwise intersections of the lines coincide with the coefficient vectors 

in (lOI) . By plugging V(^ into itself, we verify that the second iterate map equals 

(V()i)°^ = V^oV()) = -3®()i(a:, y, z) • (a: y zf . 

Hence (V(^)°^ is the identity map on all points in = 0}. Every such point 

lies in a limit cycle of length two. The points on the curve {(j) = 0} map to the 

singular points. We conclude that the ternary cubic (j) has no robust eigenpoints. 
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